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Abstract — This paper considers the problem of constructing 
a distributed direct coupling quantum observer for a closed 
linear quantum system. The proposed distributed observer 
consists of a network of quantum harmonic oscillators and it is 
shown that the distributed observer converges to a consensus in 
a time averaged sense in which each component of the observer 
estimates the specified output of the quantum plant. An example 
and simulations are included to illustrate the properties of the 
distributed observer. 

I. Introduction 

In recent years, there has been significant interest in 
controlling networks of multi-agent systems to achieve a 
consensus among the agents; e.g., see [ 1 ]—[5]. In particular, 
some authors have looked at the problem of consensus in 
distributed estimation problems; e.g., see [6], [7]. Further¬ 
more, issues of consensus have been considered in networked 
quantum systems; see [8]—[12], This work is motivated by 
the fact that it is becoming increasingly possible for quantum 
control experiments to involve the networked interconnection 
of many quantum elements and these quantum networks will 
have important applications in problems such as quantum 
communication and quantum information processing. Also, 
many macroscopic systems can be regarded as consisting of 
a large quantum network. 

In this paper, we build on the papers [13], [14] which 
considered the problem of constructing a direct coupling 
quantum observer for a given quantum system. The problem 
of constructing an observer for a linear quantum system 
has been considered in a number of recent papers; e.g, see 
[15], [16]. The control of linear quantum systems has been 
of considerable interest in recent years; e.g., see [17]—[19]. 
Such linear quantum systems commonly arise in the area of 
quantum optics; e.g., see [20], [21]. For such system models, 
an important class of control problems are coherent quantum 
feedback control problems; e.g., see [17], [18], [22]-[27], In 
these control problems, both the plant and the controller are 
quantum systems and the controller is designed to optimize 
some performance index. The coherent quantum observer 
problem can be regarded as a special case of the coherent 
quantum feedback control problem in which the objective 
of the observer is to estimate the system variables of the 
quantum plant. The papers [13], [14] considered a direct 
coupling coherent observer problem in which the observer is 
directly coupled to the plant and not coupled via a field as in 
previous papers. This leads the papers [13], [14] to consider 
a notion of time-averaged convergence for the observers. 

In this paper, we extend the results of [13] to consider 
a direct coupled distributed quantum observer which is 


constructed via the direct connection of many quantum 
harmonic oscillators. We show that this quantum network 
can be constructed so that each output of the direct coupled 
distributed quantum observer converges to the plant output 
of interest in a time averaged sense. This is a form of time 
averaged quantum consensus for the quantum networks under 
consideration. 

II. Quantum Systems 

In the distributed quantum observer problem under consid¬ 
eration, both the quantum plant and the distributed quantum 
observer are linear quantum systems; see also [17], [24], 
[28]. We will restrict attention to closed linear quantum sys¬ 
tems which do not interact with an external environment. The 
quantum mechanical behavior of a linear quantum system 
is described in terms of the system observables which are 
self-adjoint operators on an underlying infinite dimensional 
complex Hilbert space ij. The commutator of two scalar 
operators x and y on T) is defined as [x, y\ = xy — yx. Also, 
for a vector of operators x on ij, the commutator of x 
and a scalar operator y on S) is the vector of operators 
[x,y\ = xy — yx, and the commutator of x and its adjoint 
is the matrix of operators 

[: x,x *] = XX* — (x # x T ) T , 

where x# = (x^ x?j • • • x* ) T and * denotes the operator 
adjoint. 

The dynamics of the closed linear quantum systems under 
consideration are described by non-commutative differential 
equations of the form 

x(t) = Ax{t)\ x(0) = xo (1) 

where A is a real matrix in K” xn , and x{t) = 
[ Xi(t) ... x n (t) } T is a vector of system observables; 
e.g., see [17]. Here n is assumed to be an even number and 
j is the number of modes in the quantum system. 

The initial system variables x(0) = xo are assumed to 
satisfy the commutation relations 

[xj(0),x fe (0)] = 2iQ jk , j,k = l,...,n, (2) 

where 0 is a real skew-symmetric matrix with components 
Qjk- In the case of a single quantum harmonic oscillator, 
we will choose x = (xi,X 2 ) t where x± = q is the 
position operator, and X 2 = p is the momentum operator. 
The commutation relations are [q,p] = 2 i. In general, the 
matrix 0 is assumed to be of the form 

0 = diag( J, ) 


( 3 ) 


where J denotes the real skew-symmetric 2x2 matrix 


The system dynamics 0 ) are determined by the system 
Hamiltonian which is a which is a self-adjoint operator 
on the underlying Hilbert space fj. For the linear quantum 
systems under consideration, the system Hamiltonian will 
be a quadratic form H = ^x(0) T Rx(0), where R is a real 
symmetric matrix. Then, the corresponding matrix A in 0} 
is given by 

A = 2 OR. (4) 

where 0 is defined as in (J3J. e.g., see [17]. In this case, the 
system variables x(t) will satisfy the commutation relations 
at all times: 

[x(t),x(t) T ] = 2i<d for all t > 0. (5) 

That is, the system will be physically realizable ; e.g., see 

[17]. 

Remark 1: Note that that the Hamiltonian R is preserved 
in time for the system 0}. Indeed, H = ^x T Rx + ^x T Rx = 
—x T RQRx + x T RQRx = 0 since R is symmetric and 0 
is skew-symmetric. 

III. Direct Coupling Distributed Coherent 
Quantum Observers 

In our proposed direct coupling coherent quantum ob¬ 
server, the quantum plant is a single quantum harmonic 
oscillator which is a linear quantum system of the form © 
described by the non-commutative differential equation 

Zp(f) = A p Xp(t); Xp( 0) = xq p\ 

z p {t) = CpX p (t) (6) 

where z p (t) denotes the vector of system variables to be 
estimated by the observer and A p £ R 2x2 , C p £ R lx2 . It 
is assumed that this quantum plant corresponds to a plant 

Hamiltonian R p = ^x p (0) T R p x p (0). Here x p = 

where q p is the plant position operator and p p is the plant 
momentum operator. 

We now describe the linear quantum system of the form 
© which will correspond to the distributed quantum ob¬ 
server; see also [17], [24], [28]. This system is described by 
a non-commutative differential equation of the form 

Zo(^) — A 0 X 0 (t') ^ X o {0) — Xq 0 , 

Zo(t ) = c 0 x 0 (t) (7) 

where the observer output z 0 [t) is the distributed observer 
estimate vector and A p £ R" oXn °, C 0 £ R I ^ ixno . Also, 
x 0 {t) is a vector of self-adjoint non-commutative system 
variables; e.g., see [17]. We assume the distributed observer 
order n 0 is an even number with A ; ' = -rf being the 
number of elements in the distributed quantum observer. 
We also assume that the plant variables commute with the 
observer variables. The system dynamics © are determined 
by the observer system Hamiltonian which is a which is 


Qp 

Pp 


a self-adjoint operator on the underlying Hilbert space for 
the observer. For the distributed quantum observer under 
consideration, this Hamiltonian is given by a quadratic form: 
H 0 = ^a: o (0) T i? o j: o (0), where R a is a real symmetric 
matrix. Then, the corresponding matrix A a in © is given 
by 

A 0 = 2 QR 0 ( 8 ) 

where 0 is defined as in ©. Furthermore, we will assume 
that the distributed quantum observer has a chain structure 
and is coupled to the quantum plant as shown in Figure [j] 
This corresponds to an observer Hamiltonian of the form 



Z P z ol z o2 z oN 


Fig. 1: Distributed Quantum Observer. 
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where the vector of observer system variables x 0 is parti¬ 
tioned according to each element of the distributed quantum 
observer as follows 

" Zol 

Xo2 

X 0 = 

XoN 

We assume that the variables for each element of the dis¬ 
tributed quantum observer commute with the variables of all 
other elements of the distributed quantum observer; i.e., 

[x 0 i, X a j\ =0 V * 7 ^ j. 


Here, x m = 


Qoi 

Poi 


for i = 1,2,... ,N where q 0 i is the 

position operator for the zth observer element and p ol is the 
momentum operator for the ith observer element. 

In addition, we define a coupling Hamiltonian which 
defines the coupling between the quantum plant and the first 
element of the distributed quantum observer: 

H c = x p (0) T R cl x ol (0). 


Furthermore, we write 


= 


Zol 

Zo2 

ZoN 





























where 


z 0 i CoiXoi for % — 1,2, , N. 

Note that R 0 i £ R 2x2 , R C i £ R 2x2 , C 0 i £ R lx2 , and each 
matrix Rd is symmetric for i = 1, 2 ,..., N. 

The augmented quantum linear system consisting of the 
quantum plant and the distributed quantum observer is then 
a quantum system of the form ([TJ described by the total 
Hamiltonian 


Ra — Rp + Rc + Ro 

= \x a (Q) T R a x a (0) (9) 

where 



Xp 

Xol 


' R p 

Rci 





R T d 

Roi 

Rc2 

0 

X a = 

%o2 

, Ra = 

0 

R0 

Ro2 

• RcN 


%oN 





RcN RoN 


( 10 ) 

Then, using ( 0 . it follows that the augmented quantum linear 
system is described by the equations 


Xp(t) 


X p (t) 

X 0 l{t) 


Xol (t) 

Xo2{t) 

= 

Xo2 (t) 

_ X 0 N{t ) 


_ X oN (t) 


— CpX p {tf 
^o(f) - C O X O {0) 


where A a = 20R a and 
C 0 1 

C n = 


Co2 0 
0 


CoN 


( 11 ) 


We now formally define the notion of a direct coupled 
linear quantum observer. 

Definition 1: The matrices R a i, R a 2, .. R 0 n, Rci, Ra, 
..., RcN, C 0 i, C 0 2, .. C 0 at define a distributed linear 
quantum observer achieving time-averaged consensus con¬ 
vergence for the quantum linear plant 0 if the corresponding 
augmented linear quantum system O is such that 


lim — 

T-s-oo T 



1 

1 


z p (t) — z 0 {t))dt = 0. 


( 12 ) 


1 

Remark 2: Note that the above definition requires that 
the time average of each observer element output z 0 i(t) 
converges to the time average of the plant output z p (t). That 
is, an averaged consensus is reached by the observer element 
outputs. 


IV. Constructing a Direct Coupling Distributed 
Coherent Quantum Observer 

We now describe the construction of a distributed linear 
quantum observer. In this section, we assume that A p = 0 
in 0. This corresponds to R p = 0 in the plant Hamiltonian. 
It follows from 0 that the plant system variables x p (t) 
will remain fixed if the plant is not coupled to the observer. 
However, when the plant is coupled to the quantum observer 
this will no longer be the case. We will show that if the 
distributed quantum observer is suitably designed, the plant 
quantity to be estimated z p {t ) will remain fixed and the 
condition (fl2l> will be satisfied. 

We assume that the matrices R C i, R 0 i are of the form 

R ci = otifij, R oi = ujJ (13) 

where cu £ R 2 , R £ R 2 and u>i > 0 for i = 1,2 ,... ,N. 
Also, we assume that ot\ = C p . 

We will show that these assumptions imply that the 
quantity z p {t) = C p x p {t) will be constant for the augmented 
quantum system ( fm . Indeed, it follows from ( fill . (flOl) . 0 
that 

x p (t) = 2JR cl x 0 i{t) = 2Ja\/3i x 0 i(t). 

Hence, 

z p {t) = 2C p Ja\Pi Xo\{t) = 2a{ Jaiflf x 0 i{t) = 0 
since J skew-symmetric implies aj Ja\ = 0. Therefore, 

z p (t) = z p (0) = z p (14) 


for all t > 0. 

It now follows from (ITTI) that we can write 


Xol (t) 


Xol (f ) 


jPiap 

Xo2(t) 

= A 0 

Xo2{t) 

+ 2 

0 

_ X 0 n if) 


_ X oN (t) 


_ 0 


o(t)\ 


where 


A 0 = 2 


and 


= A n 


LUlJ 


X 0 l{t) 
Xo2{t ) 

X 0 iv(f) 


B 0 z p 


J 0.202 


J/320 2 L 02 J Jo3/3 3 
JP303 U13J 


J On 0n 

J/3n0n U> N J 


B n = 


2J/3i 

0 


To construct a suitable distributed quantum observer, we 
will further assume that 

ai = a = C p , 0 = —[lia and C 0 i = C p = a T (15) 


























That is, we require that 


for all i = where each fii > 0. In order to 

construct suitable values for the quantities fj n and oj t , we 
require that 

r a 


A n 


a 


B 0 \\a \\ 2 = 0. 


(16) 


m = + m+i) IHI 2 

for i = 1, 2,..., N — 1 and 

wjv = /u.jv ||a|| 2 - 


( 20 ) 

( 21 ) 


This will ensure that the quantity 


x„ = x„ - 


a 

a 


(17) 


will satisfy the non-commutative differential equation 

x 0 — A 0 x 0 - 

This, combined with the fact that 


(18) 


To show that the above candidate distributed quantum 
observer leads to the satisfaction of the condition (IT2l ). we 
first note that x 0 defined in d 1 7b will satisfy ( IT8] >. If we can 
show that 

1 f T 

lim — / x 0 (t)dt = 0 (22) 

T-s-oo T J o 

then it will follow from d 1 9b and ( IT6l ) that i fT2l i is satisfied. 
In order to establish (l22l). we first note that we can write 


where 

Ro = 


A 0 = 2 QR a 


Co; 


1 


a 

a 



UJ\I 

—fj, 2 Cta T 

~H 2 aa T 
UJ 2 I 

—fi 3 aa T 

0 


Zp 




co 31 


a 

„ T 


0 


T 

— (In olol 

—[InOIOi T 
CJjvI . 


a T 0 
0 


a 


We will now show that the symmetric matrix R 0 is positive- 
definite. 

Lemma 1: The matrix R a is positive definite. 

Proof: In order to establish this lemma, let 


' 1 ' 



1 


Xol 


Z P (19) 

Xo2 


Xo = 


1 





XoN 


will be used in establishing the condition (fl2l > for the 
distributed quantum observer. 

Now we require 


A 0 


a 

a 


= 2 


+ B 0 \\a\\ 


uJiJa — H 2 Ja||a|| 2 — fii Ja||o;|| 2 
—/r 2 Ja ||a|| 2 + u 2 Ja — /x 3 Ja||a|| 2 
— /j , 3 Ja||a|| 2 + w 3 Ja — Ja||a|| 2 

—fiiyJa\\a\\- + ojnJo. 


= 0. 

This will be satisfied if and only if 


wi -ft 2 |M| 2 — Mi||a;|| 2 
-M2IMI 2 + w 2 - M3||a|| 2 
-p 3 \\u \\ 2 +u 3 - M4||a|| 2 


be a non-zero real vector. Then 

x t 0 R 0 x 0 = u]i\\x ol \\ 2 - 2n 2 XoiCtXo 2 a + u 2 \\x 0 2 \\ 2 
- 2[i 3 x T o2 axl 3 a + to 3 \ \ x o3 \ | 2 

— 2 ^NXo N _ 1 aXg N a + at 11 cc 0 jv 11 2 
> wi||a; 0 i|| 2 - 2/r 2 ||x 0 i||||x o2 ||||a|| 2 + w 2 ||:r o2 || 2 

-2/i 3 ||a; o2 |||a; 0 3||||a|| 2 +w 3 ||a; 0 3|| 2 

— 2 ^ i n \\ x 0 n - i \\ IMIMI 2 + ccJjv ll^Coiv|| 2 

using the Cauchy-Schwarz inequality. Hence, 

x t 0 R 0 Xo > Wi||a; 0 i|| 2 - 2/i 2 ||a; 0 i||||x o2 || + w 2 ||:r o2 || 2 
—2/x 3 ||a: o2 1 ||tr 0 3 II + w 3 ||a: o3 || 2 


-fj, N \\a \\ 2 +uj n 


= 0. 


—2/iAr||a; 0 jv_i||||:E 0 jv|| + II^CoJV || 2 (23) 


where 


Pi = Pi INI 


(24) 



























for i = 1,2,..., TV. Thus. ([23} implies 


X n Ro%o 32 Ro^O 


where 


ll*ol II 

ll*o 2 || 


X n = 


and 


R„ = 


w 1 
~M 2 


0 


"M2 

W 2 

-M3 


| X 0 N | 

-M3 

w 3 


= 


-M2 M2 + M3 


-M3 


-M3 M3 + M4 


"MAT 


"MAf 

Mat 


— -Rol + -Ro 2 


where 




Mi 0 

0 0 


0 0 


> 0 


and 


Ro2 — 


M2 -M2 
-M2 M2 + M3 


-M3 
-M3 M3 + M4 


-Mat 


-Mat 

Mat 


J\T(Ro 2 ) = span{ 


}■ 


Mi M 2 


Mn 


'• —Mat 
—Mat wat 

Now the vector x 0 will be non-zero if and only if the 
vector x 0 is non-zero. Hence, the matrix R 0 will be positive- 
definite if we can show that the matrix R 0 is positive-definite. 
In order to establish this fact, we first note that ED}. (ED 
and (124} imply that 

Wi = Mi + Mi+i 
for £ = 1,2, . ..,N — 1 and 

wn = Mat- 

Hence, we can write 

Mi + M 2 M 2 


Now the matrix R a 2 is the Laplacian matrix for the 
weighted graph shown in Figure [2] 

Since this graph is a connected graph, it follows that the 
matrix R a 2 is positive-semidefinite with null space equal to 

1 


Fig. 2: Underlying weighted graph for distributed quantum 
observer. This corresponds to the observer structure shown 
in Figure [T| 


The fact that R a 1 > 0 and _R o2 > 0 implies that R a > 0. 
In order to show that R 0 > 0, suppose that x is a non-zero 
vector in Af(R a ). It follows that 

X T R 0 X = X T RolX + X T R 0 2X = 0 . 

Since R 0 1 > 0 and R a 2 > 0, x must be contained in the null 
space of R ol and the null space of /f o2 . Therefore x must 
be of the form 

‘ 1 


x = 7 


1 


1 


where 7 ^ 0. However, then 

x T R 0 \X = 7 2 Mi ^ 0 

and hence x cannot be in the null space of U„ \. Thus, we can 
conclude that the matrix R a is positive definite and hence, 
the matrix R a is positive definite. This completes the proof 
of the lemma. ■ 

We now verify that the condition (IT2} is satisfied for the 
distributed quantum observer under consideration. We recall 
from Remark |T| that the quantity \x 0 (t) T R 0 x 0 (t) remains 
constant in time for the linear system: 


= A n x„ = 2 QR n 


That is 


^x 0 (t) T R 0 x 0 (t) = ^x o (0) T R o x o (0) Vi > 0. (25) 

However, x 0 (t) = e 2eRot x o (0) and R 0 > 0. Therefore, it 
follows from ( 125} that 


V^min(Ro)\\e 20 Rot X o (O)\\ < \j\max{Ro)\\Xo 

for all 5; o (0) and t > 0. Hence, 


|e 2efl “ t || < 


^max ( Ro ) 
^min (Ro) 


( 26 ) 


for all t > 0. 

Now since 0 and R a are non-singular, 


a 20R o t 


clt = —e 


2QR 0 T n-l 























and therefore, it follows from (l26l > that 




2 QR 0 t 
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r II 2 

2T H £ 
1 

+ 2 T 


dt || 


2BR 0 T n-lci-1 




2eiioT llllfC 1 0” 1 l 


Po 1 © 


' Aroax(JR °)lli?- 1 ©- 1 l 


1 

2Ty A mini.Ro) 

+ ^ll^o 1 ©“ 1 | 


->• 0 


as T —> oo. Hence, 

lim i|| f x 0 {t)dt\\ 

T- S-OO 1 J Q 

= lim i|| f T e 2eRot x o (0)dt\\ 
T- voo J 7 0 

< lim i|| f T e 2eRot dt\\ ||x o (0)|| 
T-*-oo T J o 

= 0 . 


This implies 


1 ( T 

lim - / x 0 (t)dt = 0 
T->oo T J q 


and hence, it follows from (fTTb and ( fl9l > that 


lim — 
T—>oo T 



Also, ( IT4l) implies 


lim - / z p (t)dt 
T —>oo J Jq 


1 

1 

1 


1 

1 



Therefore, condition (fl2l > is satisfied. Thus, we have estab¬ 
lished the following theorem. 

Theorem 1: Consider a quantum plant of the form (EJi 
where A p = 0. Then the matrices li m > 0, R rl , C 0 i, 
i = 1 , 2 , ...N given as in (fl3l ). dT3T i. (l 20 t . ifTil) will define a 
distributed direct coupled quantum observer achieving time- 
averaged consensus convergence for this quantum plant. 

Remark 3: The distributed quantum observer constructed 
above is determined by choice of the positive parameters 
/ 2 i, fa, ■ ■ ■, fiN- A number of possible choices for these 
parameters could be considered. One choice is to choose 
all of these parameters to be the same as jli = u >o for 
i = 1,2,.... N where coq > 0 is a frequency parameter. 
This choice will mean that all of the oscillator frequencies 
in the distributed observer, except for the last one, will be 


the same, u>i = 2 ujq for * = 1,2,..., iV — 1 and lon = u>o- In 
order to have distinct oscillator frequencies in the distributed 
observer, we can choose jli = ioj o for * = 1,2,..., iV. This 
would yield w, = (2 i + l)wo for i = 1, 2,..., N — 1 and 
ujn = Nloq. This means that only odd harmonics of the 
fundamental frequency ujq are used. Alternatively, in order 
to obtain both odd and even harmonics of the fundamental 
frequency wo, for the case in which N is even, we can choose 

N 

M2i = = w o(— + 1 — *) >0 (27) 

for i = 1,2,..., y- This leads to 

oji = wo(IV + 1 - i) (28) 

for i = 1, 2,..., N. A similar choice can be derived for the 
case in which N is odd. 

Another possible approach is to choose the parameters 
Ai > M 2 , - - -, Ajv randomly with a uniform distribution on 

[0,c o 0 N], 

Remark 4: The proof of the above theorem relies an a 
graph theoretic argument used in the proof of Lemma Q] 
This motivates a possible extension of the result in which 
the distributed direct coupled observer corresponding to the 
weighted graph in Figure |2] is replaced by a more general 
distributed direct coupled observer corresponding to a more 
general weighted graph. 


V. Illustrative Example 


We now present some numerical simulations to illustrate 
the direct coupled distributed quantum observer described in 
the previous section. We choose the quantum plant to have 
A p = 0 and C p = [1 0]. That is, the variable to be estimated 
by the quantum observer is the position operator of the 

" 4p{t) " 


Pp(t) 


quantum plant; i.e., z p (t) = q p (t) where x p (t) = 

For the distributed quantum observer, we choose N = 5 
so that the distributed quantum observer has five elements. 
Also, as discussed in Remark 0 we choose the parameters 
Mi 5 M 2 , - - -, Miv so that jli = iujg for i = 1,2 
where loq = 1. Then the corresponding distributed quantum 
observer is defined by equations (fl3l) . (fl5l l. (f20b . d2TI> . 

The augmented plant-observer system is described by the 
equations (fTTb . ( ITOl . Then, we can write 


x a {t) = 5>(f)ir o (0) 


where 


$(f) = e Aat . 

Thus, the plant variable to be estimated z p (t) is given by 

z p (t) = ei C a $(t)x a (0) 

2N+2 

= ^2 eiC a ®i(t)x ai (0) 

i =1 


where 


C a = 


c p 0 
0 c„ 













(a) 



(b) 




Fig. 3: Coefficients defining (a) z p (t), (b) z 0 \{t), (c) z 0 2{t), 
(d) z o3 (t), (e) Zoi{t), and (f) z o5 (t). 


e\ is the first unit vector in the standard basis for R Ar+1 , 
is the ith column of the matrix <f>(f) and x a i{ 0) is 
the ith component of the vector :r a (0). We plot each of 
the quantities e 1 C a $i(t),e 1 C a $ 2 (t ),..., e 1 C a $ 2 N+ 2 (t) in 
Figure [3ja). 

From this figure, we can see that eiC a <5>i(t) = 1 and 
eiC 0 $2(f) = 0, eiC' a $ 2 (f) =0, .. e\C a $2N+2{t) = 0, 
and z p (t) will remain constant at z p ( 0) for all t > 0. 

We now consider the output variables of the distributed 
quantum observer z 0 i{t) for i = 1,2 ,,N which are given 
by 

2N+2 

Zoi (f) — ^ ^ &i-\-lC a &j(t')X a j(0') 
j=1 

where e,+i is the (i + l)th unit vector in the stan¬ 
dard basis for R w+1 . We plot each of the quan¬ 
tities ei + iCa$i(t),ei + iCa® 2 {t),...,e i+ iCa$ 2 N+ 2 {t) in 
Figures Bb) -Bf). 

Also, we can consider the spatial average obtained by 
averaging over each of the distributed observer outputs: 

^ N i N 2N+2 

z os{t ) = — 'y ' z Q i(t) = — 'y ^ 'y ] ei^.iC a ^j(t)x a j{o). 

2=1 2=1 j=1 

Then we plot each of the quantities e i+iC a &i(t), 


-ft Sill e i+lCa$2(t), ■■■ , JjYliLi e i+lC : a§2N+2(,t) in 

Figure B 



Fig. 4: Coefficients defining z os (t). 


To illustrate the time average convergence prop¬ 
erty of the quantum observer (fl2l) . we now plot the 
quantities y f 0 e i + 1 C a ®i(t)dt, y f 0 e i+1 C a ^ 2 (t)dt, ..., 
i e i+ iC a $ 2 N+ 2 (t)dt for i = 1,2,...,W in Figures 
0a)B e )- These quantities determine the averaged value of 
the ith observer output 

, ,T 2JV+2 

z oi e (T ) = - J2 e i+ iC a ^(t)x aj (0)dt 

1 Jo j=i 

for i = 1, 2,..., N. 

From these figures, we can see that for each i = 
1,2 the time average of z 0 i(t ) converges to z p (0) 
as t —> oo. That is, the distributed quantum observer reaches 
a time averaged consensus corresponding to the output of 
the quantum plant which is to be estimated. 

VI. Conclusions 

In this paper we have considered the construction of a 
distributed direct coupling observer for a closed quantum 
linear system in order to achieve a time averaged consensus 
convergence. We have also presented an illustrative example 
along with simulations to investigate the consensus behavior 
of the distributed direct coupling observer. 
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